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1. Introduction 


A general Hamiltonian system may have two different kinds of constants of motion (first 
integrals); usual ones are caused by symmetries of the Hamiltonian, while the others, 
so-called Casimir invariants, are independent of a specific choice of Hamiltonian, but are 
pertinent to the underlying Poisson algebra [DEIE]. A Casimir invariant C* is a function 
(observable) such that {C, H} = 0 for every H (here { , } denotes the Poisson bracket of 
a Hamiltonian system), i.e., C belongs to the center of the Poisson algebra. The existence 
of a nontrivial Casimir invariant {C 7^ constant) is the determining characteristics of 
a noncanonical Hamiltonian system (or, a degenerate Poisson algebra). Conversely, a 
canonical Hamiltonian system does not have a Casimir invariant; all constants of motion 
are of the first kind. 

Interestingly, we may often convert a noncanonical system into a canonical system. 
Then, what will become of Casimir invariants? Here we put the method of Clehsch 
parameterization into the perspective. The noncanonical Hamiltonian system to be 
examined is that of the ideal magnetohydrodynamics (MHD) |1]. A Hamiltonian 
formulation of ffuid/plasma mechanics is generally noncanonical when it is described 
by Eulerian variables, and has two kinds of Casimir invariants [a n El E] ; one is 
the total mass, and the others are helicities. By representing the the vector fields 
(fluid velocity and magnetic field) in terms of Clebsch parameters, we can rewrite the 
evolution equation as a canonical Hamiltonian system, and provide it with an action 
principle [H El El EDI Ell ES E31 El]- The aim of this work is to elucidate how the 
original Casimir invariants translate in the canonized system. We will show that they 
become the first-kind constants of motion corresponding to the gauge symmetries of 
the Clebsch parameterization na. The gauge transformation also conserves the action 
integral; hence the constants of motion are Noether charges, i.e., the Casimir invariants 
translate as the Noether charges corresponding to the gauge symmetry of Clebsch 
parameterization. 

We organize this paper as follows: In section]^ we start by reviewing the standard 
noncanonical Hamiltonian formulation of MHD (subsection |2.1[ ) and its canonization 
by Clebsch parameterization (subsection |2.2[ ). In section]^ we elucidate the symmetry 
of the canonized system. The Hamiltonian flow generated by the original Casimir 
invariant defines a symmetry group (subsection |3.1[ ); the symmetry turns out to be the 
gauge symmetry (redundancy) of the Clebsch parameterization (subsection |3.2[ ), which 
yields a Noether charge that reproduces the original Casimir invariant (subsection |3.3[ ). 
In subsection^, we derive an extended set of constants of motion by studying a 
generalized class of symmetries. Section concludes this work with a remark on the 
Lagrangian formalism and the relabeling symmetry. 
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2. Hamiltonian formalisms of magnetohydrodynamics 




+ -(V X B) X B, 
P 


( 1 ) 

( 2 ) 

(3) 


2.1. Ideal magnetohydrodynamics in Eulerian description 
The MHD system is described by 
d,p=-V-{Vp), 
dtV = -(V xV)xV -V 
dtB = V X {V X B), 

where p is the density, V is the velocity, B is the magnetic field, and h is the specific 
enthalpy. All variables are normalized by the Alfven unit, i.e., the energy densities 
(thermal h, kinetic pV‘^ and magnetic B‘^/2pq) are normalized by the representative 
magnetic energy density B^/po (/iq is the vacuum permeability). Here assume a 
barotropic relation h = h{p). We consider a smoothly bounded, simply connected 
domain H C We assume 

n ■ V = 0, 
n ■ B = 0, 

on the boundary dfl (n is the normal unit vector onto dfl). 

We may cast the MHD equations ([^-([^ in a Hamilton form. A general Hamilton’s 
equation is written as 

dtu = JduH, ( 6 ) 

where u E X (phase space) is a state vector, N E End(X) is a Poisson operator, and 
H E C°°{X) is a Hamiltonian. The adjoint representation of the dynamics reads 

d,F = {H,F}, (7) 

where F E C'°°(X) is an observable, and {A, 5} = {duA,J'duB) is the Poisson bracket 
induced by 77 (we denote by {u,v) the inner product of the Hilbert space X). 

For MHD, we define 

/p\ 


(4) 

(5) 


u = 


J = 


H = 


V , 



(8) 

( 0 -V- 

0 

\ 


1 

<1 

1 

1 

X 

X 

1 

X o) 

X B 

. (9) 

0 V X (o X p ^B) 

0 






(10) 


where S{p) the specific thermal energy, d{pS)/dp = h is specific enthalpy. Evidently, 
substituting ([^-(10) into (|^ reproduces the MHD equations ([^- 
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There are three Casimir invariants (hence, the MHD system is noncanonical): 


^2 = 

C'3 = 
1-1 


m 


'n 


'n 


pd?x, 

(11) 

A ■ Bd?x, 

(12) 

V ■ Bd?x, 

(13) 


where A = curP^S is the vector potential (curP^ is the inverse operator of curl, which 
is a self-adjoint operator; see [IS]); Ci is the total mass, C 2 is the magnetic helicity, and 
P 3 is the cross helicity. 


2.2. Clebsch-parameterized MHD system 

We can canonize the MHD system by parameterizing the vector helds V and B in terms 
of Clebsch parameters. Here we invoke a generalized form of Clebsch parameterization 
such that 


U = Va + bjVcj. 
i=i 


(14) 


When n is chosen to be the dimension of base space (here, n = 3), (14) gives a complete 
representation of general U by a, bj and Cj that satisfy appropriate boundary condition 
independently [IS] . From now on, the contraction rule on the indexes will be used to 
abbreviate notation. The formulation of Clebsch-parameterized MHD will be introduced 
[inidi. Here we put 


V = -V0O - —- —V/, 
P P 

B = V— X V(f)^. 

P 


(15) 

(16) 


The index ^ runs over {1, 2, 3}. 

At 9D, we assume the following boundary conditions for the Clebsch parameters: 


<1 

0 

(17) 

n ■ = 0, 

(18) 

- J 

P ’ 

(19) 

n ■ = 0, 

(20) 

n X V— = 0, 

(21) 


P 


where c^s are constant numbers, which are consistent to the boundary conditions (|^ 
and (|^ on the physical variables. Denoting 


Jr = 



( 22 ) 
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we define 


f Jc 0 0 \ 

Jc= I 0 Je 0 

0 0 X 


V 


7 


H = 


UdSx 


^p|V0o + -V7+^V07 

2 p p 


+ P^{p) + 

/ p 




Hamilton’s equation reads as a canonical system 


' <9400 = 

-W ■ V0O + - 

vy2- p-^J ■ A 

dtp = - 

-V-(pW), 


dtcP = 

-V ■ Vo^ + p-h 

J ■ v7, 

(9t7 = 

-V-(7V), 


= 

-V ■ 


. 57 ' = 

-V ■ (7V) + J 

<1 

and ^ 

J = VxB. 



We have a Lagrangian (represented by the Clebsch helds) 


L= / Cd?x = 


p(j)o + — PL ] d^x, 


eie 


by which the action is written as (denoting t = xq) 


(23) 

(24) 


(25) 


(26) 


(27) 


S= j CA‘x, (28) 

J D 

where D is a space-time domain M x hi. 

The invariants, Ci, C 2 and 6*3 are no longer “Casimir invariants” in the canonized 
system, while they must be still constants of motion. Delineation of how they are 
conserved in the canonized formalism is the aim of the next section. 


3. Gauge symmetries and the corresponding Noether charges 

3.1. Symmetry groups generated by original Casimir invariants 

Let {A,B}c = {du^A, NcducB) denotes the canonical bracket. Given a Hamiltonian H, 
a constant of motion is a functional C such that 

{G,i/}e = 0. (29) 

An inhnitesimal transformation (flow) generated by a virtual Hamiltonian G is given by 
(with small e) 


duG = eJcduJG. 


( 30 ) 
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If G and the actual Hamiltonian H commute, i.e., if G is a constant of motion, 6ug 
defines a symmetry group: (29) implies 

H{uc + Sug) - H{uc) = 0. (31) 

To simplify the notation, we will omit e in writing perturbations Suc- 

The symmetry groups generated by the “original Casimir invariants” are 

pd^x —)■ 


Gi = 


G2 = 


= 1 , 


(32) 


in 


In 


C^= V ■ Bd-^x 

Jn 

where u: = V x V. 

3.2. Gauge symmetries of the Clebsch parameterization 


500 = 

1 

<1 


6a^ = 

f V0' • B, 

(33) 

513^ = 

2V^-B, 


500 = 

(^Va' + ^V0^) •B-^V0'-n;, 


5q;^ = 

B +^V(hGu}, 

p ' p ^ ’ 


5/i^ = 

-V^-B, 

(34) 

50^ = 

-jVcfGB, 


50^ = 

-V^ ■ B + V^ -uj, 

p p ’ 



Here we show that the symmetry groups (32)-(34), generated by the original Casimir 


invariants, are gauge transformations of the Clebsch parameterization, i.e., the physical 
variables p, V and B are unchanged by the action of the group. 


Under the transformation (32), it is evident that 5p = 0 and 5B = 0, because p 


and B are independent of 0o- Also 5V = 0, because V includes 0o in the form of its 
gradient. 


Under (33), = 0 and 5B = 0 are obvious, because p and B do not depend on 

5 , and (3^. We also find that 6V vanishes: 

5(3^ 


SV = -V5(Po - 

P P 


2/ 


P 


P 


= V B ■ B 


P 


= - (V0^ • B) B 

P P P\ P ' 

= -Bx (v— X 

P VP 

= -B X B = d. 

P 


,P‘ 


1 

P 


2V— ■B\V(t)^ 
P 


( 35 ) 


Under (34), 5p = 0 is evident. 5V and 5B are calculated as follows: 


■ B - ■ n; 

p^ p^ J p^ 


6V = - V 










Gauge Symmetries and Noether Charges in Glebsch-Parameterized MagnetohydrodynamicsJ 


1 

P 


- ( Vj ■B]Va^ 


+ —V ( -Va^ ■ B - -0^ ■ a; 
P \P P 

P\P P 

+ —V ■ B 


P 


P 


= - (v— ■ b\ Va^ - - (Va^ ■ B) V— 
P\P ) P P 


- (V/ ■ Lj) V— - - ( V— ■ LJ I V0 


p 




1 

p p 
£ 


P 




- ( V— ■ S ) - - (V/ • B) V— 

p V p / p p 

= -S X X V— + V(/)^ X V—^ + -a; X 

P V P Pj P 

= -(Bxu) + ujx B) = 0, 

P 

P P 


V— X V(/)' 
P 


(36) 


= - V X 


- ■ S) 

p V p / p p 




= V X -B X V— X V0' 


P 


P 


= V X 


-S X B 1 = 0. 

P 


(37) 


S'.,?. Noether charges 

To show that Ci, C 2 and 6*3 are Noether charges of the canonized system, we have yet 
to prove that the gauge transformations (the symmetry groups) conserve the action. 

The Lagrangian density C consists of two parts; the canonical 1-form and Pi. 
Because the transformations (32)-(34) are the gauge transformations, it is obvious 
that PL is unchanged. On the other hand, the canonical 1-form may be changed by 
a gauge transformation. However, the change is an exact differential, thus the action 
integral is changed only by a constant. In fact, an arbitrary Hamiltonian flow duc 
(generated by a virtual Hamiltonian G) conserves the canonical 1-form. Let us denote 
duc/e = itc = deUc = Jcdu^G. For pi = (p,/r,/3) and g* = (0 o,Q!, 0), we observe 


Pi9' = Gipiq') + hcJc 'lie. 


(38) 


and 
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= -2ucdu,G+ de{ucdu,G), 


(39) 


which implies that is the exact differential. Explicitly, the action of (33) yields 

5C = p6^o + 


iP ■ 

h' At-, 


= -dt{A-B) + W ■[ + A - V ^ 


P 


1 / . 

P 11 


P 


X A 


(40) 


thus 5C is an exact differential. Similarly, by the action of (34), we obtain 
dC = p6^o + P^da^ + dp^a^ + + 5 / 3 ^^ 

= -dt {V ■ B) 




v|dy 

p 


X V’ 


(41) 


P P 

where V' = V + V0o = -p^/pWa^ - /^VpV/. 

Let us briefly review how a Noether charge is derived by an action-invariant 
transformation. We consider a first-order Lagrangian such that L = £(g\ di,q^)d^x, 

which yields the Euler-Lagrange equation dgi£ — d^, [dg^giC) = 0. Suppose that L has a 
symmetry: dC = (with some vector A'^) under a transformation g* —)■ g* -|- 5g*. For 
a solution q’‘(t,x) of the Euler-Lagrange equation, we obtain 

dC = dq'^^q^C d{dyq'')dQ^C 

= dq^ {^g^C - dy {da^qiC)) + {dq^dg^giC) 

= {dq^dg^grC) . ( 42 ) 

For a perturbation dq'’ of a symmetry group, we may write 

d,{dq^dg^g.C~N^)=d,r = ( 43 ) 

where B is a Noether current, and I^d^x = [dq^dgiC — A°) d^x is the corresponding 
Noether charge. 

Applying this procedure to the canonized MHD Lagrangian and the gauge 
symmetries, we obtain 


d(j)o = 1 —)■ / pd'^x = Cl, 

Jn 


S- V 0 ■ (jO 


d(j)o = 

1 

<! 

to 

r 


da^ = 


B 

^ A- 

■ Bd^x = 

djd^ = 

2V^ ■ 

p 

B 

Jn 


d(f>o = 

a 

l> 

■ B - 

1 

<1 

doA = 

—-VA 
p , 

'■B + 

4 • uj 


dp^ = 

-V^ . 

p 

B 



dcj)^ = 


B 



djd^ = 


B + \ 

p 



(44) 

(45) 


-)■ 


V' ■ Bd^x. 


(46) 
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The third Noether charge (46) is equal to the cross helicity C 3 ; using the boundary 
condition of B, 


V ■ Bd?x 


in 


V-Bd^x= / V(l>o-Bd^x 
Jn 

= / V ■ (00^) 

Jn 

= n ■ (poBd^x = 0. 
Jan 


(47) 


Remark 1 (Gauge transformation changing the action) A transformation con¬ 
serving a Lagrangian maps a solution of the eguation of motion to another solu¬ 
tion. However, the converse is not necessarily true In fact, there is a gauge trans¬ 
formation of the Clebsch parameterization which does not conserve the action. Let 
6(j)o = —• B, 5a^ = • B, and 6f3^ = 0oVy ■ B. These perturbations yield 

a gauge transformation: 

5V = -V(50o - 

P P 


= ^ {V(p^ ■B)V^- — (V^-B)V(p 


P 

= —B X R = 0, 
P 


P P 


7P 


P 


(48) 


It is also obvious that B and p are unchanged. However, we obtain 

5C = p6^o + p^Sd^ + 5 / 3 ^^ 


= (Po(-dt{A-B) + V- ( ^^^B 


A - V ( 


X A 


, (49) 


which is not an exact differential. 


We have elucidated a beautiful relation, mediated by the constants Ci, C 2 and G 3 , 
between the Eulerian noncanonical formalism and the Clebsch-parameterized canonical 
formalism; in the former, the constants are the Casimir invariants that foliate the phase 
space, and in the latter, they are the Noether charges that identify symmetry groups 
reflecting the redundancy of the Clebsch parameterization (see Fig.[^. 


3 . 4 . Generalized constants of motion 

The Clebsch parameterization enables us to discribe local conservation of mass and 
helicity. We define two sets of functions Gi and G 2 such that 

Gi = {a\dta + V ■ Va = 0} , (50) 

G 2 = {Al^tAT V-(AV) =0}. (51) 

From (26), (p^ G Gi and p, fy G G 2 . We easily verify 

A, p G G 2 , s, f G M —y sX T fp G G 2 , (52) 
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gauge symmetry 


Figure 1. The relation between Casimir invariants and Noether charges. Left: In the 
noncanonical formulation in terms of Eulerian physical variables, Casimir invariants 
foliates the phase space. Orbits are constrained on manifolds that are intersections of 
the energy contours (level-sets of H) and the Casimir leaves (level-sets of C). Right: 
By the Clebsch parameterization, the system is embedded in a larger phase space with 
redundant degree of freedom; the Hamiltonian group action generated by the Noether 
charges delineates the symmetry built in the Clebsch parameterization. 


a, 6 G Gi —)■ /(a, b) G Gi for any smooth function f, (53) 

A,r] G G2 —)■ — G Gi, (54) 

V 

a G Gi, A G G 2 -)■ aA G G 2 , (55) 

a, 6, c G Gi — )■ Va ■ (V6 x Vc) G G 2 . (56) 


When X E G 2 , Ad^x is a constant of motion. From (54), = yf /p G Gi and using 

(53) and (55), pf{cr^, G G 2 ; hence we hnd a new constant of motion pficr^, 0^)d^x 
(/ is an arbitrary smooth function). From (55) and (j5^,cr*V0-^ ■ (Va^ x Va^) G G 2 . 
Using ( [ 5 ^ , and summing up the terms such that i = j and k = I, we obtain A B E G 2 . 
From ( [^ , /(cr^, (f)^)A ■ B E G 2 , thus /(cr^, <f^)A ■ Bd^x is a constant of motion. 

Since /(cx^, 0^) is a scalar function that is constant along the streamlines of the 
flow V, we may regard it as a weighting function co-moving with the plasma. Choosing 
/(cr^, 0^) to be a support function on a co-moving volume element, p/(<x^, 0^)d^x and 
In (j)^)A ■ Bd^x measure local mass and helicity contained in the volume element. 
These generalized constants of motion are related to the following gauge symmetries: 




in 


S<t'o = f- f,/f, 
fa' = /f, 

= -Pfi, 


(57) 
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f S^„ = -l^fiA ■ B - l^V^‘ ■ (/S + V X (/A)), 
^ J Sa‘^ iflA ■ B + iV4>‘ ■ (/S + V X (/A)). 

[ SP’ = -/'A . B + V ■ (l^ (/B + V X (/A))) , 


(58) 


where /f denotes the derivative of / with respect to a^, and the one with respect to 
(j)^. For example, 6 V under the transformation (57) is 


SV = -V (^/ - + /'V/ - dv/j' 

= -V/ + f^v4>' + fly- 

P 

= -V/ + V/ = 0, (59) 


and 6 p = t) and 6 B 


0 are obvious, thus (57) 


is a gauge transformation. 


4. Conclusion and remarks 

We have shown, in the context of MHD, that Casimir invariants of a noncanonical 
Hamiltonian system are converted, upon canonization by an alternative parameteriza¬ 
tion (here, Clebsch parameterization) of the state vector, into Noether charges pertinent 
to the gauge freedom (redundancy) of the new parameterization. While the conclusion 
delineates a simple relation between different formalisms of dynamical systems, the 
specihc forms of the symmetries, embodying the concrete relation, are rather compli¬ 
cated. Since the Clebsch parameterization implies a nonlinear variable transformation, 
the symmetries are not apparent (excepting some simple ones; cf. [15] )• Emphasizing 
another aspect of the problem, we may say that we could unearth, guided by the known 
Casimir invariants, the general gauge symmetries of the Clebsch parameterization. 

There are some different methods of canonization. A well-known method that 
applies for various fluid-mechanical systems is the usage of Lagrangian variables to 
represent the dynamics of fluid elements; the Lagrangian variables label the initial 
position of each fluid element. Formulating an action by a Lagrangian of Newtonian- 
mechanical continuum, we obtain a canonical systems of Hamilton’s equation of 
motion ; see [T71 [18] for the original formulation of the Lagrangian action principle 
of MHD. In the canonized system of Lagrangian variables, the Casimir invariants are 
translated differently. The total mass Ci and the magnetic helicity C 2 are now trivially 
conserved, because they are the integrals of the attributes of fluid elements; the local 
density and magnetic helicity are not dynamical variables, but are constants bound 
to each fluid element (like charges of particles in particle systems). Only the cross 
helicity C 3 is a nontrivial hrst integral; it is now the Noether charge pertaining to the 
relabeling symmetry that represents the arbitrariness of providing fluid elements with 
Lagrangian labels [I211I91120]. Since Clebsch parameters may be viewed as the Eulerian 
counterparts of Lagrangian labels (or, the pull-back along the Cauchy characteristics 
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of the fluid velocity V ; see m), it is natural that the gauge symmetry of the Clebsch 
parameters (present canonization) and the relabeling symmetry of the Lagrangian labels 
(Lagrangian canonization) yield the same constant of motion. 
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